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Abstract. We study the evolution and emission of circumbinary disks around close classical T Tauri binary systems. High 
resolution numerical hydrodynamical simulations are employed to model a system consisting of a central eccentric binary star 
within an irradiated accretion disk. A detailed energy balance including viscous heating, radiative cooling and irradiation from 
the central star is applied to calculate accurately the emitted spectral energy distribution. 

Numerical simulations using two different methods, the previously developed Dual-Grid technique with a finite difference 
discretization, and the Smoothed Particle Hydrodynamics method are employed to compare the hydrodynamical features and 
strengthen our conclusions. 

Physical parameters of the setup are chosen to model the close systems of DQ Tau and AK Sco. Using the self-consistent 
models, we are able to fit the observed spectral energy distributions by constraining parameters such as disk mass, density 
profile and radial extension for those systems. We find that the incorporation of irradiation effects is necessary to obtain correct 
disk temperatures. 

Key words, accretion disks - binaries: spectroscopic - hydrodynamics - methods: Numerical 



1. Introduction 

From observations it is known that a large fraction of young 
stars are in binaries. The formation of binary systems by close 
int eractions in a star forming region was studied numerically 
by bate et al.l J2002I) . Typically, a circumbinary accretion disk 
surrounding the double star is formed during this process. Here 
\ we concentrate on a later phase and model the hydrodynamic 
evolution of the circumbinary disk around close binary stars 
having a separation of only a fraction of an AU. This is inspired 
by observations of the spectroscopic binaries DQ Tau and AK 
Sco for which there is also a spectral energy distribution (SED) 
av ailable. The latter sy stem AK Sco was reanalyzed recently 

bv lAlencar et alJ d2003l) . 

In a previous paper fciiinther & Klevll2002l) we modeled 
the evolution of circumbinary disks and compared fully de- 
veloped circumbinary disks and their properties with ob- 
ser vational data. That work exten ded previous computations 
bv lArtvmowicz & Lubowl Jl996ll and iRozvczka & Laughlinl 
ill 9971) by solving explicitly a time dependent energy equation 
including the effects of viscous heating and radiative cooling. 
We investigated both the structure and dynamics of the disk as 
well as the gas flow in the close vicinity of the binary star. To 
that purpose we utilized a newly developed method which en- 
ables us to cover the whole spatial domain called the Dual-Grid 
method. For the first time we performed long-time integration 
of the complete system covering several hundred orbital peri- 



ods of the binary and compared the properties of the evolved 
systems with observational data such as spectral energy distri- 
butions in the infrared and optical bands and accretion rates 
estimated from observed luminosities. 

In this paper we extend the model by accounting for irra- 
diation effects which are important for passive accretion disks 
such as the disks of DQ Tau and AK Sco, according to observa- 
tiona l data. Additionally, r efined observational data is available 
from lAlencar et all d2003l) for the AK Sco system which allows 
to better constrain the parameters of the system by our numer- 
ical simulations. 

In the next section we present the equations that are be- 
ing solved and the refined radiative balance m odel. The overall 
layout of the physical model is the same as in lGunther & Klevl 
d2002l) . Then we proceed with some remarks on the generation 
of the spectra from the simulation data (Sect.[3)- After that in 
Sect.0]we describe the two-dimensional simulation setup and 
proceed to the new results in Sect. [5] We summarize and con- 
clude in Sect. [6] 

2. Equations 

The evolution of the disk is given by the two-dimensional (r, tp) 
evolutionary Navier-Stokes equations for the density X, the ve- 
locity field u = (u r , Up), and the temperature T. In a coordinate- 
free representation the equations read 



Send offprint requests to: R. Giinther, 

e-mail: rguenth@tat . physik . uni - tuebingen . de 



— + V • (Eh) = 



(1) 



2 



Richard Gunther, Christoph Schafer & Wilhelm Kley: Evolution of irradiated Circumbinary Disks 



ffLu 
~dt 

<9Ee 
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+ V • [(Se + P)u] = -Im ■ <X> - V • (F - u ■ T) 



(2) 



(3) 



Here e = c v T + 1 /2m 2 is the specific total energy, P is the ver- 
tically integrated (two-dimensional) pressure P - RYT '//u with 
the midplane temperature T and the mean molecular weight jj., 
which can be obtained by solving the Saha rate equations for 
Hydrogen dissociation and ionization and Helium ionization. 
F and T are the radiative flux and viscous stress tensor, respec- 
tively. 

The gravitational potential <S> generated by the binary stars 
is given by 

GMj 
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for \r -n\ > R* 



for \r - r,| < R t 



(4) 



where r, are the radius vectors to the two stars, and R* is the 
stellar radius assumed to be identical for both stars. The second 
case in Eq. gives the potential inside a star with radius 
having a homogeneous (constant) density. 

2.1. Viscosity and Disk Height 

The effects of viscosity are contained in the viscous stress ten- 
sor T. Here we assume that the accretion disk may be de- 
scribed as a viscous medium driven by some internal turbu- 
lence which we approximate with a Reynolds ansatz for the 
stress tensor. The components of T in differ ent coordinate s ys- 
tems are spelled out explicitly for example in lTassoull(ll978l) .A 
useful form for disk calcula tions conside ring angular momen- 
tum conservation is given in lKlevI dl999h . 

For the kinematic she ar viscosity v we use an o--model 
JShakura & Sunvae\ll973l) of the form 



v = ac s H , 



where c s and H are the local sound speed and vertical height, 
respectively. The local disk height H(r) is computed from the 
vertical hydrostatic equilibrium which yields 



H(r) 



Z 
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(6) 



This can be split into the single star disk heights given by 

Hi{r) - c s J ^qm, ■ The sound speed is given here by c s = 
RT/fi. 

2.2. Radiative balance 

The influence of radiative and viscous effects on the disk tem- 
perature is treated in a local vertical balance. The balance equa- 
tion for the internal and radiative energy considering these ef- 
fects reads 



d(Lc w T) 
dt 



- Gdiss - Giad + 6s 



(7) 



where <2diss and g la d denote the viscous dissipation and radia- 
tive losses. Qstars is the irradiation heating by the stars, which is 
considered for the optically thick regions of the circumbinary 
disk. This inclusion of irradiation effects is done in a simple 
geometric fashion due to the lack of explicit treatment of the 
vertical structure of the disk in these two-dimensional calcula- 
tions. For the viscous dissipation we use the vertically averaged 
expression 
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where T is the viscous stress tensor. For the radiative transport 
we use 

dF z 
dz 

where Fo is the flux vector in the z — plane. Here we consider 
only the losses in the vertical direction, i.e. Fo = 0, a stan- 
dard approximation in accretion disk theory. Integration over 
the vertical direction yields 

grad = 2F rad = 2(T B T 4 eff (10) 

with the local effective (surface) temperature r e ff. Following 
iHubenvH 1990ll this temperature is related to the midplane tem- 
perature T by the LTE gray solution for the vertical structure of 
accretion disks. Simplifying this relation we get 
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The optical depth is calculated via dRuden&Pollackll99lh 



T = — 

2 



(12) 



using an interpolated opacity x(p, T) adapted from 
Lin & Papaloizouldl985l) . 



2.2.1 . Irradiation 



(c\ For the heating from the stars we use 



(13) 



where T* is the effective blackbody temperature of the stars, R, 
is the stellar radius and r, are the distances of the stars to the 
local patch. The factor A honors the effective surface the flux 
operates on. This is obtained by reconstructing a surface from 
the radially varying pressure scale height. Self-shadowing is 
taken into account by setting A to zero for back-faces. 

Numerical tests have shown that with the inclusion of irra- 
diation effects we observe oscillatory behavior of the tempera- 
ture near to the edge of the gap du e to self-shadowing effects 
as also found bv lDullemondlj2000l) . But we can overcome this 
specific instability by using appropriate time sub-stepping for 
the radiative balance. This suggests that if the radiative balance 
is imposed using a time relaxation process the self-shadowing 
instability does not occur. This may be either because the cool- 
ing timescale is not within the region of linear instability or 
because of non-linear effects exposed by the relaxation mecha- 
nism. 
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Fig. 1. Visible setup of the AK Sco disk system using the esti- 
mated inclination of 63°. The disk shape is modeled according 
to pressure scale height, gray scale coding is logarithmic sur- 
face density in g/cm 2 . The disk extent is 12 AU. 

3. Spectra Generation 

Spectral energy distributions from the disk emission can be 
computed by integrating the flux of the blackbody radiation 
from the disk surface. This is done decoupled fro m the numeri- 
cal sim ulation as a postprocessing step. We follow lAdams et all 
Jl988l) . and compute the flux at the observer via 

cos i r Rd r^" 

F v = J J o B v [T(r, cp)} (l - e~ T ™) rdtpdr (14) 

where i is the inclination of the disk, D the distance to the ob- 
server and t the line-of-sight optical depth through the disk, 
which can be approximated using the opacity k, r(r, tp) = -f^f-p ■ 
The surface disk temperature T is estimated using the optical 
depth of the disk as in the radiative balance process. 

This integration is limited to regions outside of the stel- 
lar cores as temperatures and densities inside the cores are 
not modeled correctly. The stellar emission is accounted for 
by adding two blackbody spectra for appropriate effective stel- 
lar temperatures. To compare with observational data all of 
the generated spectra are reddened according to an extinction 
A v suggested by the cited references using the method from 
ICardelhet al.1dl989h . 

4. Simulation Setup 

We start with an axisymmetric circumbinary disk around the 
binary system with an initial gap of a width that is determined 
by simulations including the gap formation process. This re- 
sults in gap widths about three to four t imes the binary separa- 
tion ra dius consistent with results from lArtvmowicz & Lubowl 
(11994 . This initial configuration is then evolved for several or- 
bital periods, so the circumstellar environment can form and 
the whole system settles into a quasi-periodic state. 

Initial density and temperature distributions f or the disk 
are taken from power-law fits of iMathieu et alJ J 19971) and 




Fig. 2. Linear surface density (g/cm 2 ) of the inner circumbi- 
nary disk of the DQ Tau system after 50 orbital periods. Scaling 
is in AU. The spiral structure in the circumbinary disk induced 
by the binary is clearly seen. 



lAlencar et al1(l2003l) . This initial surface temperature is relaxed 
to a midplane temperature satisfying Eq. (II It . 

The stellar temperatures and radii specify the flux observed, 
and that transfered into the disk. These parameters cannot be 
fitted independently but need to be fixed for the irradiation 
of the disk. This is why we use stellar radii as suggested by 
IMathieu et ail (1 19971) and lAlencar et alJd2003l) and fit the effec- 
tive temperatures T t according to the observed spectral energy 
distributions by assuming simple blackbody spectra. Direct 
matching of the stellar spectra is a good approximation as can 
be seen in Fig.^ where we display the AK Sco disk with the bi- 
nary stars projected using the estimated inclination of 63°. This 
demonstrates that occultation of the stars by the circumbinary 
disk is not to be expected. 

The overall azimuthally averaged circumbinary disk shape 
does not change significantly on timescales of the orbital period 
of the binary system. Induced by the gravitational torques of 
the central binary a nice a spiral structure develops in the inner 
regions of the disk as can be seen in Fig. [2] 

Also one can observe periodic formation of spiral arms 
from the edge of the circumbinary disk going down to the cir- 
cumstellar environment (Fig.|3|l connected to the binary orbital 
period. The quasi-periodic state of the circumstellar environ- 
ment of the eccentric system is characterized by the following 
phases: 

1 . Going from periastron to apastron, material is accreted to 
the circumstellar environment from the circumbinary disk 
and the remaining material in the gap. 

2. Around apastron phase circumstellar envelopes/disks have 
been formed. 



4 



Richard Gunther, Christoph Schafer & Wilhelm Kley: Evolution of irradiated Circumbinary Disks 




Fig. 3. The spiral arms onto the circumstellar environment for 
the DQ Tau system after 50 orbital periods. The gray scale 
follows logarithmic surface density annotated with iso-lines. 
Plotted is the region inside the circumbinary disk gap which is 
located at 0.4 AU. The binary is in apastron phase. 

3. Approaching periastron, these circumstellar disks are torn 
off due to gravitational torques exerted on the disks by the 
stars, and a part of the mass is accreted onto the stars. 

This leads to perio dic accretion event s which is both confirmed 
from observations dBasri et all 19971) showing periodic bright- 
ening events and from numerical simulations. A more detailed 
discussion on the disk d ynamics and accretion can be found in 
iGiinther & Kievl(l2002l) . 

4.1. Dual-Grid and SPH 

We additionally have performed Smoothed Particle 
Hydrodynamics (SPH) test calculations which cover only 
the inner part of the simulation domain to corroborate the 
results of the Dual-Grid technique. The initial setup consists of 
two circumstellar disks in Keplerian motion around the indi- 
vidual stars in apastron. Similar to the grid based calculations 
we evolve this system for a couple of orbital periods to reach 
a quasi-stationar y state. For deta iled descrip t ions o f the SPH 
method see, e.g.. lBena Jl99d) and lMonaghai] ( ll992t) . The SPH 
simulations include only radiative cooling and viscous heating 
and not the irradiation from the stars, as these calculations 
are primarily intended to support the results obtained by the 
Dual-Grid technique. 

5. Results 

Grid based calculations have been performed on two different 
regions of the whole system. First, high resolution (1 10 x 204) 
calculations extending only over the circumbinary disk (0.4— 
12 AU for AK Sco and 0.3-80 AU for DQ Tau) have been per- 
formed to evolve these regions in higher resolution and for a 
long time. Second, high resolution ( 1 80 x 247 for the r — <p grid, 
101 x 101 for the Cartesian grid) calculations extending over 




Fig. 4. Circumstellar disks logarithmic density structure for 
Dual-Grid simulation (left) and SPH simulation (right) while 
approaching periastron. The binary separation is about 0.1 AU. 



the circumstellar disks and up to the very inner region of the 
circumbinary disk (up to 1 .2 AU for AK Sco and up to 1 .0 AU 
for DQ Tau) have been run to allow comparisons with the SPH 
calculations. Finally, low resolution (128 x91, 37 x37) simula- 
tions of the whole system using the Dual-Grid technique have 
been performed to show that indeed these regions decouple in 
the timescales we are interested in. 

5.1. Dual-Grid and SPH 

The Dual-Grid technique was developed to overcome limita- 
tions of using cylindrical coordinates for accretion disks, where 
the whole domain incl uding the origin needs to be covered 
JGiinther & Klevll2002t) . To demonstrate that the interpolation 
scheme used between the Cartesian and the cylindrical grid is 
accurate enough for modeling close binary circumstellar disks, 
we performed comparison runs using SPH. 

Simulations are set up to initially contain two circumstellar 
disks around the individual stars with a mass of 10 -8 M and a 
disk profile according to £ oc r~ 3 ^ 2 . This system is then evolved 
for a few orbital periods to truncate the disks to the right size 
and go into a quasi-stationary equilibrium. Dual-Grid calcula- 
tions were performed on a 180 x 247 sized r — <p grid and a 
101 x 101 sized Cartesian grid, while the SPH calculations in- 
strumented 150 000 particles. 

Fig- |H shows the structure of the disks for Dual-Grid and 
SPH calculations after the same evolution time. In both sim- 
ulations one can see similar circumstellar material features, 
namely 

- spiral structure inside the circumstellar material, 

- spiral waves going off the circumstellar region, 

- a bar-like connection between both circumstellar regions. 

Thus we can conclude that with the Dual-Grid calculations we 
see the same features across the origin as with the SPH cal- 
culations. The Dual-Grid technique is applicable to simulation 
of circumbinary systems including the origin and the inner cir- 
cumstellar regions. 

5.2. Irradiation Effects 

For motivating the inclusion of irradiation effects, we first an- 
alyze the individual contribution of the different terms in the 
radiative balance (Eq. [7}, and the pressure work pV ■ u, to the 
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Table 1. Parameters for DQ Tau jMathieu et al J 19971) . 



Fig. 5. Contribution of the different heating/cooling terms to 
the rate of temperature change (DQ Tau after 50 orbital peri- 
ods). Only the inner part of the system is displayed. The irradi- 
ation irr and the radiative cooling -cool match exactly beyond 
0.4 AU. 
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Fig. 6. Contribution of the different heating/cooling terms to 
the rate of temperature change (DQ Tau after 50 orbital peri- 
ods). Plotted is the outer circumbinary disk part of the system. 
The irradiation by the central stars and the radiative cooling 
match exactly (topmost dashed line). 



temperature change rate. The azimuthally averaged rate of the 
temperature change at quasi-stationary equilibrium of the DQ 
Tau system is plotted in Fig.|5]for the inner regions of the disk 
and Fig. [6] for the outer regions. Here, vise denotes the rate 
from the viscous heating, irr the rate from the irradiation pro- 
cess, cool the emission contribution and pdV the rate from the 
pressure work. Note, that the inner edge of the disk is approxi- 
mately at r — 0.4 AU for this model. 

As one can easily infer from Fig.|6]the circumbinary disk is 
dominated by irradiation effects (which in fact match the cool- 
ing rates exactly), and as such qualifies as a passive disk. The 
effects due to viscous heating and pressure work can be ne- 
glected here. 
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Table 2. Parameters for AK Sco JAlencar et all2003h . 



In the inner optically thin part the pressure work is domi- 
nating due to the dynamic motion in this region caused by the 
gravitational torques exerted by the eccentric binary. One can 
also see that the emitted flux in this region is negligible com- 
pared to the flux emitted by the circumbinary disk. 

5.3. Circumbinary Disks 

We model the two close spectroscopic binary systems with 
circumbinary disks, DQ Tau and AK Sco, which are both of 
T Tauri t ype. The phys i cal pa rameters for th e systems are 
taken f rom lAlencar et alJ d2003l) for AK Sco and lMathieu et alJ 
ll 19971) for DQ Tau. The relevant information for the simula- 
tions has been summarized in Tables^andlJ] Both stars are as- 
sumed to have the same effective temperature T* and the same 
radius R„. 

For both the DQ Tau and the AK Sco system we shall see 
missing flux at around 5 x 10 13 , resp. 10 14 Hz which can be 
accounted to the modeling of the T Tauri type stars emission as 
a simple blackbody. 

5.3.1. DQ Tau 

From full simulations including the gap formation process we 
infer an initial gap with a radius of 0.4 AU for the DQ Tau 
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Fig. 7. Comput ed SEP of DQ Tau. Cr osses show observational 
data taken from lMathieu et al.ldl997l) . 
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Fig. 8. Azimuthally averaged surface temperature of the equi- 
librium DQ Tau circumbinary disk. 



system. The stars combined effective blackbody temperature 
is fitted manually to 3500 K to match the observations. Starting 
from the parameters prese nted in Tableland power-law fits 
from M athieu et alJ i 19971) we adjust the initial density distri- 
bution power-law index, the disk mass and its extension to best 
match the observed SED after going to quasi-stationary state. 

The SED of DQ Tau can be fitted best with a thin disk ex- 
tending to 80 AU which has an unusually fiat surface density 
distribution £ oc r -°- 25 and a total disk mass of 7 x 10~ 4 M Q . 
This is less mass a nd a larger disk radius as obtained by 
iMathieu et alJ dl997h . After fifty orbital periods, the equilib- 
rium surface temperature distribution of the circumbinary disk 
follows T(r) = 155 K r~ 42 as can be seen in Fig. [8] 

5.3.2. AK Sco 

In Fig. 16 of lAlencar et alJ J2003I) the authors propose a cir- 
cumbinary disk mass of Md = 0.02 M and a disk extend of 
Rd = 12 AU for the best fit to the observed SED. The surface 



Fig. 9. Computed SED of AK Sco with paramet ers and obser- 
vation al data (crosses with error-bars) taken from lAlencar et al] 
J2003I) . 
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Fig. 10. Computed SED of the AK Sco system with best-fit 
parameters . Crosses with error -bars show observational data 
taken from lAlencar et al.l d2003l) . 



density distributio n is that of a minimum mas s solar nebula as 
in the used model (Chiang & Goldreichl 19971) . so E(r) ~ r _1 . 
In Fig.[9]the emission of a model evolved with these parameters 
imposed initially is displayed. 

Clearly visible is the excess flux of the model in the in- 
frared (around 5 x 10 13 Hz) which we can reduce by lowering 
the disk mass and slightly flattening the disk density profile. 
Using a disk mass of Md = 0.01 M and an initial density pro- 
file according to E(r) ~ r~ L0 we obtain a spectrum as seen in 
Fig. [101 

For the AK Sco system the stars combined effective black- 
body temperature is best fitted to the observations by assuming 
T» = 6500 K. In a quasi-stationary state after about fifty orbital 
periods of the binary the surface temperature can be fitted to 
the power-law T{r) = 264 K r~ 0A1 for r > 2 AU and a surpris- 
ingly steep linear behavior in the region between one and two 
AU(r(r) = 617 K-215 K r) as shown inFig.lTTl This is due to 
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Fig. 11. Azimuthally averaged surface temperature of the equi- 
librium AK Sco circumbinary disk. 



the different opacities in the dense region of the disk and com- 
pared to the case of DQ Tau three orders of magnitude higher 
maximum surface density which is related to the AK Sco disk 
being much smaller and more massive. 

We note that our models are not able to fit the higher flux 
observed in the SED of both systems in the region between the 
disk and stellar contributio ns. This m ay be attributed to line 
emission features dAlencar et alj|2003l) which cannot be mod- 
eled by our simulations. 

5.4. Circumstellar Disks 

For the region inside the circumbinary disk gap we see a tiny 
amount of warm gas and streamers feeding circumstellar disk 
like structures from the inner circumbinary disk edge. Due to 
the low mass this material does not contribute to the continuum 
part of the observed spectral energy distributions but rather 
would show up in lines and the UV part of the spectra. 

Also, the observed periodic brightening in the light-curve 
from the systems originate here, which is believed to orig- 
inate from periodic accretion events. These periodic events 
of acc re tion have been confirmed by lArtvmowicz & Lubow] 
Jl996j) . iRozvczka & LaughlirJ d 19971) and iGiinfher & Klevl 
J2002I) through numerical simulations. 

The circumstellar disk like structures around DQ Tau form 
through accumulation of material in the stellar gravitational po- 
tential. Without removing any material through an accretion 
process, disk profiles according to Fig. I12lform. They cannot 
be identified with a classical accretion disk due to shape and 
size as one can see from Fig.ll3lwhich shows the circumstel- 
lar material and velocity distribution. They rather would form 
sort of an envelope, but this remains to be investigated in three- 
dimensional calculations. 

6. Conclusion 

We model the circumbinary disks of DQ Tau and AK Sco and 
their circumstellar environment including emitted spectra by 



Fig. 12. Profiles of the circumstellar material along the con- 
necting line of the binary at apastron in the system of DQ Tau 
(stellar radii 0.008 AU). 




Fig. 13. Circumstellar material in the DQ Tau system after 41 
orbital periods, corresponding to Fig.^] Gray scale coding is 
velocity magnitude annotated with vector glyphs, isolines are 
equally spaced logarithmic surface density. The stellar cores 
are the white circles. 



applying two-dimensional hydrodynamical simulations of high 
accuracy and resolution. The models include a self-consistent 
approximative vertical energy transport including irradiation 
from the stars using detailed opacities and equation of state. 
The numerical resolution we have reached by now is so fine 
that the individual stars in close binaries are already resolved 
by typically about 10x10 grid-cells for high-resolution calcu- 
lations. 

We find that with suitable variations of 

1 . the initial surface density profile 2 ~ r~ d , 

2. the circumbinary disk mass Mi 

3. and the circumbinary disk extent 

the spectral energy distributions of DQ Tau and AK Sco can be 
fitted well to the observational data only if irradiation effects 
are included. 
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The circumbinary disk settles into a quasi-stationary state Ruden, S. P. & Pollack, J. B. 1991, ApJ, 375, 740 
after around 50 orbital periods of the binary with a temper- Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 
ature distribution following approximately T(r) oc r -° 45 for Tassoul, J. 1978, Theory of rotating stars (Princeton Series in 
larger radii in both systems. We find that irradiation plays the Astrophysics, Princeton: University Press, 1978) 
dominant role for the heating balance in these outer regions. In 
the inner part of the disk just beyond the inner gap, the addi- 
tional heating due to the tidal effects is important, particularly 
for small disks around tight binaries, such as AK Sco. 

For the present model parameter we do not find any in- 
dication for a previously observed self- shadowing instabi lity 
of the disk due to the stellar irradiation dDullemondll2000l) . A 
proper numerical treatment of the strongly non-linear cooling 
is required to prevent artificial numerical instabilities, which 
is achieved by using a time-relaxation scheme as opposed to 
directly solving for the equilibrium. 

Future models of this kind need to include a more detailed 
treatment of the stellar boundary layers as well as a more de- 
tailed model of the accretion process. We believe that doing full 
radiative transfer in two dimensions will not improve our re- 
sults, but going to three dimensions with full radiative transfer 
included will possibly enable us to obtain more accurate spec- 
tra from the circumbinary disk gap region with its thin material. 
We then may be able to compare in more detail with observed 
light curves and phase dependent spectra. 
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